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Abstract 

This preliminary report is the foreword of a series dedicated to stochastic deforma- 
tions of curves. Problems are set in terms of exclusion processes, the ultimate goal 
being to derive hydrodynamic limits for these systems after proper scaHngs. Here, 
only the basic ASEP system on the torus Z/NZ is analyzed. The usual sequence 
of empirical measures, converges in probability to a deterministic measure, which 
is the unique weak solution of a Cauchy problem. The method presents some new 
features, letting hope for extensions to higher dimension. It relies on the analy- 
sis of specific partial differential equations involving variational calculus. Namely, 
the variables are the values of functions at given points and their number becomes 
infinite. 



Keywords Exclusion process, hydrodynamic limit, martingale, Cauchy prob- 
lem, weak solution, distribution, parabolic operator. 

1 Preliminaries 

Interplay between discrete and continuous description is a recurrent question 
in statistical physics, which in some cases can be answered quite rigorously 
via probabilistic methods. In the context of reaction-diffusion systems, this 
is tantamount to studying fluid or hydrodynamics limits. Number of ap- 
proaches have been proposed, in particular in the framework of exclusion 
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2 Model definition 



processes, see e.g. [25, [Z| ISE], ^Hj and references therein. As far as fluid 
limits are at stake, all these methods have in common to be limited to sys- 
tems for which the stationary states are given in closed product forms, as 
far as hydrodynamic limits are concerned, or at least for which the invariant 
measure for finite (the size of the system) is explicitly known. For instance, 
ASEP with open boundary can be described in terms of matrix product form 
(a sort of a non-commutative product form) and the continuous limits can 
be understood by means of brownian bridges (see (H]). We propose to adress 
this question from the following different point of view: starting from dis- 
crete sample paths subjected to stochastic deformations, the ultimate goal 
is to understand the nature of the limit curves when increases to infinity. 
How do these curves evolve with time, and which limiting process do they 
represent as t goes to infinity (equilibrium curves)? Following ^31 and [T^ . 
we will try to give some partial answers to these questions in a series of 
papers. 

This first study is mainly dedicated to the ASEP model. The mathematical 
approch relies on the analysis of specific partial differential equations involv- 
ing variational calculus. A usual sequence of empirical measures is shown 
to converge in probability to a deterministic measure, which is the unique 
weak solution of a Cauchy problem. Here variables are the values of some 
function at given points and their number becomes infinite. In our opinion, 
the method presents some new features, which let hope for extensions to 
higher dimension. 

A future concern will be to establish a complete hierarchy of systems of 
hydrodynamic equations , whose steady state will help to describe non-Gibbs 
states. 

2 Model definition 

2.1 A stochastic clock model 

The systems we will consider can typically be describded as an oriented path 
embedded in a bidimensional manifold: A^ steps of equal size, each one being 
chosen among a discrete set of n possible orientations drawn from the set 
{2k-K/n,k = 0...n — 1} of angles with some given origin. The stochastic 
dynamics which is applied consists in displacing one single point at a time 
without breaking the path, while keeping all links within the set of admissible 
orientations. In this operation two links are displaced. This constrains quite 
strongly the possible dynamical rules, which are given in terms of reactions 
between consecutive links. 



2.2 Examples 
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We have 

2.2 Examples 

(1) The simple exclusion process 

The first elementary and most sudied example is the simple exclusion pro- 
cess: this model, after mapping particles onto links, corresponds to a one- 
dimensional fluctuating interface. Here we have a binary alphabet and letting 
= r and X"^ = r, the set of reactions simply rewrites 

A- 

rf ^ f r, 
A+ 

where are the transition rates for the jump of a particle to the right or 
to the left. 

(2) The triangular lattice and the ABC model 

Here the evolution of the random walk is restricted to the triangular lattice. 
Each link (or step) of the walk is either 1, e^*'^/^ or e^*'^/^, and quite naturally 
will be said to be of type A, B or C, respectively. This corresponds to the 
so-called ABC model, since there is a coding by means of a 3-letter alphabet. 
The set of transitions (or reactions) is given by 

AB ^ BA, BC h CB, CA ^ AC, (2.1) 

p+ q+ r+ 

where there is a priori no symmetry, but we will impose periodic boundary 
conditions on the sample paths. This model was first introduced in in 
the context of particles with exclusion, and for some cases corresponding to 
the reversibility of the process, a Gibbs form for the invariant measure was 
given in |12| 

3 Hydrodynamics for a one-dimensional asymmetric exclusion 
process [ASEP] 

As mentioned above, we aim at obtaining hydrodynamic equations for a class 
of exclusion models. The method, although relying on classical powerful 
tools (martingales, relative compactness of measures, functional analysis), 
has some new features which should hopefully prove fruitful in other contexts. 
The essence of the approach is in fact contained in the analysis of the popular 



4 3 Hydrodynamics for a one-dimensional asymmetric exclusion process [ASEP] 



ASEP model, presented below. We note the difficulty to find in the existing 
literature a complete study encompassing various special cases (symmetry, 
total asymmetry, etc). Some proofs will only be sketched, and the related 
results presented as claims or even conjectures. 

Consider N sites labelled from 1 to A^, forming a discrete closed curve in the 
plane, so that the numbering of sites is implicitly taken modulo A^, i.e. on 
the discrete torus G*^' = Z/A^Z. In higher dimension, say on the lattice Z'^, 
the related set of sites would be drawn on the torus {JLjNTL)^ . 

We gather below some notational material valid throughout this section. 

• 71 stands for the real line. C'^[0, 1] is the collection of all real- valued, k- 
continuously differentiable functions defined on the interval [0, 1] , and 
M. is the space of all finite positive measures on the torus G [0, 1). 

Co°(i^) is the space of infinitely differentiable functions with compact 
support included in K. 

• For S an arbitrary metric space, 7^(S) is the set of probability measures 
l?s[0,r] is the space of right continuous functions z : [0, oo] — > S with 
left limits and t ^ Zt- 

• For i = 1, . . . , A^, let A^^\t) and B^-^\t) be binary random variables 
representing respectively a particule or a hole at site i, so that, owing 
to the exclusion constraint, Af\t) + B^'^^t) = 1, for all 1 < i < A^. 
Thus {A('^)(t) = {Af\t), . . .,A'^^\t)),t > 0} is a Markov process. 

• ri'^' will denote the generator of the Markov process A'^'(t), and 
^j^' = (t(A'^'(s),s < t) is the associated natural filtration. 

• Our purpose is to analyze the sequence of empirical random measures 



when A^ — > oo, after a convenient scaling of the parameters of the 
generator Q'-'^K The probability distribution associated with the path 
of the Markov process G [Oi^]i for some fixed T, is simply 

denoted by Q^'^l 

As usual, one can embed G'^^ in G, so that a point i € G'^' corresponds 
to the point i/N in G. Hence, in view of H3.H1 . it is quite natural to let 
the sequence Q'^' be defined on a unique space T>m[0,T], which becomes a 
polish space (i.e. complete and separable) via the usual Skorokod topology, 
as soon as ^A is itself Polish (see e.g. ^HI, chapter 4). Without further 
comment, ^A is assumed to be endowed with the vague product topology, as 
a consequence of the famous Banach-Alaoglo and Tychonoff theorems (see 




(3.1) 



e.g. EBlIIEl). 
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Let (t)a,4>b be two arbitrary functions in C^[0,1] and define the real- valued 
positive measure 



^t''^'[0a,0fe] = exp 



(3.2) 



which is a functional of (paj'pb- For the sake of brevity, the explicit depen- 
dence of A^'^\t), B!:'^\t), Zj:'^^ [4>a, (j'b] On N, t, (f>, will be omitted wherever the 
meaning remains clear from the context: for instance, we often shall simply 
write Ai,Bi or Z^'^K Also Z'-'^'' stands for the process {Z^'^\ t > 0}. 

A standard powerful method to prove the convergence (in a sense to be 
specified later) of the sequence of probability measures introduced in 
consists first in showing its relative compactness, and then in verifying the 
coincidence of all possible limit points (see e.g. ^HI- Moreover here it suffices 
to prove these two properties for the sequence of projected measures defined 
on V-jilOjT] and corresponding to the processes {Zl'^^[(f)a, <Pb]ii > 0}, since 
the functions belong to C^[0, 1]. 

Let us now introduce quantities which, as far as scaling is concerned, are 
crucial in order to obtain meaningful hydrodynamic equations. 

^ ' 2 ' (3.3) 

^fliN) = Xab{N)-Xta{N), 

where the dependence of the rates on N is explicitly mentioned. 

Theorem 3.1. Let system hS. have a given asymptotic expansion of the 
form 

' X{N) = XN'' + o{N^), 

(3.4) 

^I2{N) =^,N + o{N), 

where X and n are fixed positive constants. [As for the scaling assumption 
Jg.^l ), the random measure log Zj.'^'' is a functional of the underlying Markov 
process, in which the time has been speeded up by a factor N"^ and the space 
shrunk by N^^J. Assume moreover the sequence of initial empirical measures 
log^Q^'', taken at time t = 0, converges in probability to some deterministic 
measure with a given density p{x,0), so that 

/■I 

AN) 



lim log.^Q = / [p{x,0)(f>a{x) + {1 — p{x,0))<f>b{x)]dx, in probability, 

N^oo Jq 

(3.5) 

for any pair of functions (j)a,4'b G where K £ TZ is a compact con- 

taining the interval [0, 1]. 
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Then, for every t > 0, the sequence of random measures nl converges in 
probabilitys, as N ^ oo, to a deterministic measure having a density {p{x, t) 
with respect to the Lebesgue measure, which is the unique weak solution of 
the Cauchy problem 



/o Jo 



de{x,t) 



dx 



dxdt 



[p{x, T)9{x, T) - p{x, 0)9{x, 0)] dx, 

(3.6) 

where holds for any function d G C^([0, 1] x [0,T]). 

//, moreover, one assumes the existence of ^-^p^, then i.V. 6)) reduces to a 
classical Burger's equation 

dp{x,t) aVOM) , r, „ . 

— A — h ^[1 - 2p{x,t)\- 



dt 



dx"^ 



dx 



Proof. The sketch of the proof is spread over three main subsections, re- 
ferred to hereafter as PI, P2 and P3. 



PI [Existence of limit points: sequential compactness] As usual in prob- 
lems dealing with convergence of sequences of probability measures, our very 
starting point will be to establish the weak relative compactness of the set 
{logZ^^^A^ > 1}. Some of the probabilistic arguments employed in this 
paragraph are in a way classical and can be found in good books, e.g. |25|ll9|. 
although for simpler models. 

The process 

^(AT) def ^(AT) _ ^(iV) _ f ^(N)^Z^N)^^S (3.7) 

^0 

is a bounded {.Fj^'j-martingale. Using the exponential form of Zl'^^ together 
with classical stochastic calculus (see e.g. jl^, chap. 3, page 93), it follows 
that 

[yW] ^=1:' ([/W)2 _ ^*(^J7W[(Zf ))2] - 2Zf' ])ds (3.8) 
is also a bounded real martingale. 

Prom A'^Pit) + Bf\t) = 1,V1 < i < A^, on sees that Z^""^ is mainly a 
functional of the sole function if^xy = <Px ~ 4'y = ~'4'yx, up to a constant 
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uniformly bounded in A^. Hence, setting 



def , fi + l 



tpxy 



Ky{i,N) = X,y{N) 

we have 



exp( ^AV'..,^- 



xy = ab or ba, 



where 



(3.9) 
(3.10) 



By using the exclusion property, a straightforward calculation in equation 
(|3.1()|l allows to rewrite H3.8|l in the form 



(3.11) 



where the process i?^^' is stricly positive and given by 



<> 



Ki>{N) 



The integral term in I|3.1H1 is nothing else but the increasing process associ- 
ated with Doob's decomposition of the submartingale (C/^'^')^. 



The folllowing estimates are crucial. 
Lemma 3.2. 



0{1), 



(3.12) 
(3.13) 



Proof. We will derive H3.12|l by estimating the right-hand side member of 
equation H3.inp . 

Clearly, Aip^yl^j^^ = -^''Pxy(^'k^ + ^(iT^)' where ip' denotes the derivative 
of tp. Then, taking a second order expansion of the exponential function and 
using definitions (|3.3p and (EHJ, we can rewrite (j3.1()p as 



LI 



(JV) 



N 
X{N) 



AiAi^ 



(3.14) 
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The first sum in H3.14|l is uniformly bounded by a constant depending on ^jJ. 
Indeed, < 1 and ip € C^[0, 1], so that ip' is of bounded variation. 

As for the second sum coming in H3.14|l . we have 



Then the discrete Laplacian 

admits of the simple form 

A*4i±i) - A^^i) = (3.15) 

where ^p" denotes the second derivative of ip. 

By (IS3I), A(iV) = XN^ + o{N'^), so that iTTra implies 

(3.16) 

which concludes the proof of (|3.12|l . The computation of R['^^ leading to 
(|3.13p can be obtained via similar arguments. ■ 

To show the relative compactness of the family Z'-^\ which here, by separa- 
bility and completeness of the underlying spaces, is equivalent to tightness, 
we proceed as in pHI by means of the following useful criterion. 



Proposition 3.3 (Aldous's tightness criterion, see [4j). A sequence {X'-^^} 
of random elements ofVulOjT] is tight (i.e. the distributions of the {X'-^^} 
are tight) if the two following conditions hold: 

(^) 

lim limsupP[||X<'^'|| > a] = 0, (3.17) 



N 



where = supjXf | 



(a) For each e,r], there exist a 6o and an Nq, such that, if 5 < 6q and 
N > Nq, and if T is an arbitrary stopping time with r + 6 < T , then 

P[\Xl^^,-Xr\>e] <v. (3.18) 



Note that condition hS. 1 is always necessary for tightness. 
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We will now apply Lemma, [T2l to equations (j8.7|l and (|3.1ip . the role of X'-^^ 
in Proposition 13.31 being played by Z*^'. 

Observe that, by the uniform boundedness of Zj-^\ condition (|3.17p is im- 
mediately verified. 

To check condition H3.1Sjl . rewrite (|3.7|l as 

- zr = ui:i - ur + f^^\^-\zi-^ds. (3.19) 



The integral term in (|3.19|l is bounded in modulus by K5 [where -fC is a 
constant uniformly bounded in N and V'] and hence satisfies (|3.18p . We are 
left with the analysis of C/j^''. But, from (|3.1ip . H3.13p and Doob's inequality 
for sub-martingales, we have 



E[(^ 



(AT) 

t+5 



p 



sup I U^ 

t<T 



(iV)| 



{iv)\21 



> e 



< 



t+5 



{Z'f^fR[''^ds 



C 



< 



4C_ 



(3.20) 



(iV) 



almost 



where C is a positive constant depending only on ip. Thus 
surely as ^ oo. This last property together with assumption H3.5p yield 
(|3.18p and the announced (weak) relative compactness of the sequence Z^^' . 
Hence, the sequence of probability measures Q^^', defined on T>m[0,T] and 
corresponding to the process /ij'^', is also relatively compact: this is a con- 
sequence of classical projection theorems (see for instance Theorem 16.27 in 
[TH]). We are now in a position to state a further important property. 

Let Q the limit point of some arbitrary subsequence Q^^''\ as Uk — > oo, 
and Zt = lim„j,^oo z|"''\ Then the support of Q is a set of sample paths 
absolutely continuous with respect to the Lebesgue measure. Indeed, the 
application fit — > sup^<2^ log Zt is continuous and we have the immediate 
bound 



sup log Zt 

t<T 



< 



[\Mx)\ + IM^Wx, 



which holds for all ipa, tpb ^ C^[0, 1]. Hence, by weak convergence, any limit 
point Zt has the form 



Zt[(t)a,<pb] = e-xp / [p{x,t)(t)a{x) + {1 - p{x,t)4>b{x)]dx 



(3.21) 



where p{x, t) denotes the limit density (a priori random) of the sequence of 
empirical measures z/^™'"'* introduced in H3.ip . 
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P2 [A functional integral operator to characterize limit points] This 
is somehow the Gordian knot of the problem. Relying on the above weak 
compactness property, our next result shows that any arbitrary limit point 
Q is concentrated on a set of trajectories which are weak solutions of a 
functional integral equation (IFE). 

First, by IjH.Tp . (pHUl and we obtain at once 



dt 



,2 ylN) 



(3.22) 



It is worth remarking that H,S.22|1 should be written, strictly speaking, as 
a stochastic differential equation, which is well-defined since indeed all the 
underlying probability spaces emanate from a families of interacting Poisson 
processes. 

Replacing for a while the quantities (f'aij^) and 4>b{jj) by variables xf^^ and 



(iV) 



respectively, (13.2211 becomes 



(JV) 



dt 



d Z 



3 7(N) 



i+l 

(3.23) 



where we have put 



exp 



exp 



(AT) , (JV) 



(JV) 



N 



(JV) WV, I 



(JV) 



We shall rewrite (|3.23|1 in the operator form 

-^=cr[zri (3.24) 

remarking in the present setting that, for each finite N, acts on the 
function space C^[— |(/'|]^'^, where 

|<A| = sup {\Mz)\,\Mz)\), (3.25) 
2e[o,i] 

7(N) 



and p is an arbitrary positive number, as Zj: is analytic with respect to 

.(j^ 
-t 



{(l)a{-),4>b{-)}- The operator is of parabolic type, but then in the wide 
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sense, since here one can check the quadratic form usually associated with 
the second order derivative terms is non definite, see e.g. [9J). 

The key point will be to show that any limit point Zt ^= lim z'f'^^ satisfies 

an IFE, obtained by studying the second order linear partial differential 
operators along the sequence Uk ^ oo. 

To carry out the analysis of the limit sum coming in (|,S.2,Sp (which is a priori 
intricate), we propose a general approch, which aims at proving first that Zt 
is a weak solution (or distributional in the sense of Schwartz) of a Cauchy 
type operator. The line of argument will be sketched below. 

Beforehand, for the sake of shortness, it will be convenient to define the 
following cylinder sets, for p = 1, 2 . . ., 

uf'='[-\cp\,\cp\rx[o,t], up'='[-\<p\,mp. 

Introduce the operator O/^^ which is the adjoint of in the Lagrange 
sense, so that, for every function h G C^{Uf''), 

~(^) def ^ , ^^2 V- d"^ [a^yji , N)h] d'^ [ayg^ji , N)h] 

4|f |^ + ^(iv)[;,], (3 26) 

Definition 3.4 (see e.g. [24j, Part 1, section III. 5). A function g ^ L2 is 
said to be a weak (or distributional) solution of the Cauchy problem C^^^g = 
if, for all h G C^iU^"), 

gC[''^[h]dudt = 0, (3.27) 

where in the integral u denotes an arbitrary point in U'^^ . 

Multiplying equation H3.23|l by an arbitrary function h G Cq^{U.^'^), for fixed 
T arbitrary positive, and then integrating twice by parts, we obtain, in 
agreement with (j3.27p . 



U2N 



dudt = 

(3.28) 

[(^w _ uP)h{u,T) - Z^^^/i(n,0)]dn. 



A brute force analysis of the adjoint operator could lead to a dead-end. A 
preliminary step will be to exploit carefully the estimates obtained in Lemma 
13.21 This is the content of the next lemma. 
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Lemma 3.5. The following partial differential equation holds 

d{zr-un ^ .... fi 



(3.29) 



where the term lis i?J modulus uniformly bounded by ^, C being a 

constant depending only on and tp" . 

Proof. Immediate from equations (|3.14|l and HH.1()|1 . ■ 

Starting from Lemma f,S.5| we will present the two global guidelines of a func- 
tional approach, called Gl and G2. Basically, it relies on partial differential 
equations, whose variables are functions taken at points of the torus. We 
think this might well extend to larger dimensions, although this assertion 
could certainly be debated. 

Gl Intermediate reduction to an almost sure convergence context. This 
can be achieved by means of the extended Skohorod coupling theorem 
(see Corollary 6.12 in jTHj), which in brief says that, if a sequence of real 
random variables (^fc) is such that limfc^oo /A:(Cfc) = /(O converges in 
distribution, then there exist a probability space V and a new random 
sequence ^fc, such that = & and limfc_oo /fc(Cfc) = fiO, almost 
surely in V, with C = Here this theorem will be applied to the 
family Zj^"'''^ which thus gives rise a new sequence denoted by y^^"'*-' 
in the sequel. This step is in no way obligatory, but just a matter of 
taste. Indeed, one could still keep on with weak convergence context 
and use Alexandrov's portmanteau theorem (see e.g. pHI) whenever 
needed. 

G2 For each finite A^, we can consider the quantities 

as constant parameters, while the x^^' 's are free variables. This is 
clearly feasible, choosing for instance (pai-), (t>b{-) in the class of polyno- 
mials of degree at least 3A^. Also, from now on, the functions 4>a and 
(ph will be supposed to belong to C^{K), for some compact K £ TZ 
containing the interval [0, 1]. 
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Then, according to Gl, we rewrite (|8.29p as 



dt 



(3.30) 



where A\- is viewed as an operator of parabolic type with constant coefH- 
cients and domain C^iplj, 

def dg 



Aria] 



dt 



I f dg ^ dg 



N 



d'g 



remembering that ipab = 4>a — fpb- The term C>(^^^ iii l|3.3r)|l stands for an 
operator having a neghgible range for ^ oo. 

Let denote the adjoint of ^J'^^ Then 



dh 



N 



If dh dh \ d'^h 



2 V dxl"^^ ^^i+i 



dx\'''>dx[l{ 



i+1 



(3.31) 



and, for any /i'^' € C^{UP), we have 



(N) 



(3.32) 



where 5u in (|3.32|l represents the differential volume element 

5u = dxfUxf\..dx';:^' = 5</,„ (1)^0,(1) . . . <5</>,(l). 



The next step is to make a suitable choice of the function h in (|3.32p in 
order to extract a meaningful information on the limit operator, as ^ oo. 
Keeping in mind that the random variables Y^'^^ are defined on the implicit 
probability space V introduced in Gl, we state the following result. 
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Lemma 3.6. [Claiml] Let Yt = lim a.s., and let W G C^{K) be a 

fixed compact space. Then, for a class of test functions k properly chosen in 
a subset ofC^{W x [0,r]), we have 

T 



(3.33) 



dt J Ytm]m{-)m-) = 

[k{^{.),T)YTm]-H^i-),0)Yom]]5H-), 



where 



dk{cp{.),t) 



dt 



dx. 



Moreover IIS. .V/j)) yields the Cauchy problem posed in i.V. 

Proof. We only sketch the main hnes of argument. 

• The almost sure convergence of Y^^^ and C/j^' , respectively to Yt and 
[by iHOl ]. will yield (TOil provided that in iTO^ the functions /i'^' 
are properly chosen, to ensure the existence of the limit sums coming 
in iTOTTl . as ^ oo. 



AN) 



Setting X — ^u.^ ,^^2 ) 
that it suffices to take in H3.32|l 



(^(N) (iV) 



, x^'), he reader can convince himself 



h^-^[^,t)=kix'-\t), 



where k € Cq^(W x [0, T]), in which case convergent Riemann sums 
are obtained in I|3.3H1 . 

From Skohorod's coupling theorem, Yf does satisfy an equation of the 
form H3.2H1 . Hence, we can write the following functional derivatives 
(which are plainly of a Radon-Nykodym nature) 

BYt 



d^{.) - 

19<A2(.) 



p{;t)Yt, 
= p\;t)Yt. 



To derive (EH, one has to pick out k from a class of convolution test 
functions, depending on some parameter e and properly converging in 
the space of Schwartz distributions. 
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P3 [Uniqueness] The problem of uniqueness of weak solutions of the 
Cauchy problem (|3.6|l for nonlinear parabolic quation is in fact already solved 
in the literature. We refer the reader for instance to f9j for a wide bibliogra- 
phy on the subject. 

To conclude the analysis of Theorem l3.11 it suffices to switch back to Z^'^^ [(f>a, 4>b] , 
which converges in distribution to Zt[(j)a,4>b]- Hence, the random measure 
/xj^' converges in distribution to a deterministic measure, which is a peculiar 
situation impling also convergence in probability. ■ 



4 The n-species model 

We will state a conjecture about hydrodynamic equations for the n-species 
model, briefly introduced in section m\ in the so-called equidiffusion case, 
precisely defined hereafter. 

The n-species system is said to be equidiffusive whenever there exists a con- 
stant D, such that, for all pairs {k,l), 



Letting 



def , XkliN) 

aki = hm log ■ 



we assert that the following hydrodynamic system holds [Claim2]. 



dpk 
dt 



D 



dx"^ dx 



l^k 



, fc = 1, . . . , n. 



The idea is to apply the functional approach presented in this paper. 
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